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INTRODUCTION 

The  Luzon  Strait  is  a  key  place  in  the  East  Asian  climate  change  study.  The  climate  condition  of  the 
South  China  Sea  is  significantly  controlled  by  the  Luzon  Strait  transport  (LST)  (Qu  et  ah,  2004;  Liu  et 
ah,  2008).  Variability  in  the  Kuroshio  origin  around  the  strait  is  mostly  caused  by  the  LST,  and 
significantly  affects  the  oceanic  variability  around  Japan  and  Korea  in  the  downstream.  In  spite  of  the 
importance  of  LST  measurement,  most  of  the  previous  observations  are  based  on  hydrography  (Nitani, 
1972;  Shaw,  1991;  Chen  and  Huang,  1996;  Yaremchuk  et  ah,  2009),  and  the  direct  measurement  of 
current  is  still  limited  in  the  observational  period  and  spatial  coverage  (Y.  Yuan  et  ah,  2008;  Liang  et 
ah,  2008). 

The  ocean  acoustic  tomography  (OAT)  is  a  powerful  tool  that  can  estimate  the  structure  of  sound 
speed  (mainly  proportional  to  temperature)  and  current  in  the  ocean  (Munk  et  ah,  1995).  While  a  lot  of 
effort  has  been  devoted  to  sound  speed  structure  measurement,  application  of  OAT  to  current  structure 
measurement  is  still  restricted  to  the  coastal  seas  (Park  and  Kaneko,  2000;  Yamaguchi  et  ah,  2005  and 
Lin  et  ah,  2005).  The  OAT  has  never  been  applied  in  oceans  characterized  with  strong  internal  tides 
and  waves  like  the  Luzon  Strait  (Niwa  and  Hibiya,  2004).  This  may  be  caused  by  that  the  sound 
transmission  in  such  oceans  interacts  strongly  with  these  tides  and  waves  and  is  unexpectedly  disturbed 
not  to  obtain  an  accurate  result  of  temperature  and  current. 

A  long-term  OAT  experiment  for  current  measurement  is  attempted  at  the  northern  part  of  the  Luzon 
Strait  where  the  Kuroshio-related  current  and  the  strong  internal  tides  and  waves  occur.  Four-month 
observation  results  are  presented  together  with  the  new  clock  correction  method  and  compared  with  the 
moored  upward-looking  acoustic  Doppler  current  profder  (ADCP)  data. 

OBSERVATION 

A  middle-range  OAT  experiment  composed  of  three  mooring  stations  (Tl,  T2  and  T3)  was  carried  out 
with  range  about  36km  at  the  northern  part  of  Luzon  Strait  during  April  to  September  2008  (Fig.  1). 

The  floor  depth  ranges  from  900m  to  1700m  around  the  tomography  site.  The  800Hz  organ  pipe 
transducer  (Engineering  Acoustic  Inc.)  is  attached  at  depth  800m  on  the  subsurface  mooring  line  for  all 
stations  and  transmits  one 
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Fig.  1  Location  map  of  the  observation  site.  The  bathymetric  chart  is  shown  in  a  magnified  scale  on 
the  right  side  of  the  figure.  Tl,  T2  and  T3  are  the  mooring  stations  equipped  with  the  OAT  system. 
T1  and  T2  also  equip  the  upward-looking  moored  ADCP  at  the  top  of  the  mooring  line. 


period  (about  two  minutes)  of  the  12*  order  M  sequence  every  eight  hours.  The  eight-hour  data  are 
smoothed  through  the  14-day  running  mean  to  filter  out  the  tidal  variation  and  focus  on  stable  ocean 
currents  like  the  Kuroshio. 

The  transducer  is  a  narrow  band  one  with  a  frequency  range  of  25Hz  while  the  conventional  acoustic 
tomography  experiments  always  have  used  broad-band  transducers.  This  forces  us  to  shift  to  the  cycle- 
per-digit  value  of  24  from  4  suitable  for  the  conventional  experiment.  The  power  of  sound  transmission 
can  be  considerably  reduced  in  the  operation  of  narrow-band  transducer.  This  provides  a  big  advantage 
for  the  long-term  operation  of  the  subsurface  tomography  system  by  reducing  the  battery  consumption. 
The  deployment  and  recovery  works  were  performed  onboard  the  R.V.  Dong  Fang  Hong  of  the 
Ocean  University  of  China.  The  four-month  data  were  successfully  acquired  for  the  stations  Tl  and  T2 
while  the  T3  mooring  could  not  be  recovered  by  the  troubles  which  may  be  caused  by  fishing  nets. 

Note  that  the  top  float  of  the  mooring  line  is  located  at  depth  150m  for  T3  in  contrast  to  Tl  and  T2  at 
depths  600m  and  300m,  respectively. 

The  moored  upward-looking  ADCP  (RDl-SC  Narrow-Band  for  Tl  and  RDI-SC  Workhorse  Long 
Ranger  for  T2)  is  attached  at  the  top  of  the  mooring  lines  Tl  and  T2  for  comparing  with  the  OAT  data. 
The  working  frequencies  of  ADCP  are  150kHz  and  75kHz  for  Tl  and  T2,  respectively.  Unfortunately 
the  Tl-ADCP  was  suddenly  stopped  in  one  and  half  months  after  the  deployment  because  of  the 
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system  trouble.  The  raw  ADCP  veloeities,  aequired  with  bin  lengths  of  8m  and  20m  for  T1  and  T2, 
respeetively  are  also  smoothed  by  the  two-day  running  mean  (Tl)  or  40-h  Butterworth  low-pass  filter 
(T2)  to  remove  the  semi-diurnal  and  diurnal  tides.  A  number  of  CTD  casts  were  done  around  the 
tomography  site  at  the  deployment  and  recovery  cruises  to  know  the  sound  speed  fields. 

Clock  correction 

Special  attention  is  here  paid  on  the  accurate  measurement  of  ocean  current  which  requires  highly 
precise  clock.  The  OAT  system  is  equipped  with  a  timing  module,  composed  of  GPS  (global 
positioning  system)  clock,  quartz  clock  (SeaScan,  Time  Base  Module)  and  Rubidium  clock  (TEMEX 
TIME,  MCERS-01).  The  quartz  clock  is  always  operated  as  an  internal  clock  with  calendar  during  the 
experiment.  The  internal  clock  is  synchronized  with  the  GPS  clock  immediately  before  the  deployment 
into  water  and  the  drift  of  the  internal  clock  is  measured  in  comparison  with  the  GPS  clock 
immediately  after  the  recovery  from  water.  The  Rubidium  clock  is  operated  during  the  first  day  after 
the  deployment  and  compared  with  the  internal  clock  to  estimate  the  drift  rate  per  second.  The  internal 
clock  is  corrected  with  this  drift  rate  before  the  next  transmission.  This  procedure  of  clock  correction  is 
well  operated  during  the  whole  observation.  However  the  resulting  clock  drifts  make  several  hundred 
milliseconds  much  greater  than  the  expected  ones  of  ten  milliseconds.  This  clock  drift  is  too  big  to 
estimate  ocean  currents  in  a  36km  range  with  a  satisfactory  accuracy  of  about  5cm/s.  The  new  clock 
correction  method  is  here  proposed  to  replace  the  timing  module  method.  When  an  acoustic  transducer 
is  located  near  the  axis  of  the  underwater  sound  channel  (EiSC),  some  acoustic  rays  advance  near  the 
use  axis  with  no  current.  This  requires  that  travel  time  differences  are  zero  for  these  rays.  By  using 
this  assumption,  the  time-wise  growth  of  clock  correction  (At)  is  formulated  with  the  following 
regression  curve; 

A^  =  -0.001 +  0.0028f"  -  0.000 if'  +  0.0102f '  +  0.1960f  +  0.2553 
t  ={t-  Mean{t))/  STD{t) 

where  t  denotes  the  time  in  unit  of  hour  elapsed  from  the  deployment  of  the  system  into  water  and  At 
the  clock  correction  in  unit  of  second  at  time  t.  The  Mean(t)  and  STD(t)  denote  the  mean  and  standard 
deviation  of  t,  respectively  The  data  scatter  around  the  above  line  with  a  STD  of  10.3ms,  which 
corresponds  to  the  velocity  of  4. 1  cm/s.  This  velocity  provides  an  error  bar  for  the  accuracy  of  the  new 
clock  correction  method,  proposed  here. 

Ray  simulation  and  the  sound  transmission  data 

The  sound  transmission  process  between  the  stations  Tl  and  T2  is  simulated  by  the  ray-tracing 
method,  in  which  only  the  reflection  and  refraction  of  sound  are  considered.  This  is  performed  by 
using  the  range-averaged  sound  speed  profde,  estimated  from  all  the  acquired  CTD  data.  The  ray  paths 
with  bottom  reflection  numbers  less  than  two  are  divided  into  three  groups,  shown  with  the  ray 
diagram  in  Pig.  2.  The  travel  times  for  the  rays  are  plotted  against  the  launch  angles  at  lower  middle  of 
the  figure.  The  first  arrival  group  (Ray-1)  makes  travel  times  around  23.9s.  Also  the  arrival  times  for 
the  second  (Ray-2)  and  third  (Ray-3)  arrival  groups  are  around  24.0s  and  24.1s,  respectively.  Ray- 
passes  the  depth  range  of  5  00- 13  00m.  Ray-2  and  Ray-3  reach  the  200m  layer  at  the  northern  two-third 
of  the  section  while  they  are  quite  different  at  the  southern  one-third  of  the  section  where  Ray-2 
draws  a  downward  sloping  line  and  Ray-3  an  upward  convex  curve. 


3 


The  transmission  signals  at  the  stations  T1  and  T2  are  phase-modulated  with  the  12*^  order  M  sequence 
of  different  codes.  The  received  data  at  the  stations  T1  and  T2  are  cross-correlated  with  the  same  codes 
of  the  12*'^  order  M  sequence,  used  in  the  transmission  at  T1  and  T2  to  increase  the  signal-to-noise  ratio 
(SNR).  The  travel  time  of  acoustic  ray  between  the  stations  T1  and  T2  can  be  determined  from  the 
significant  arrival  peaks  which  exist  in  the  correlation  waveform.  As  shown  typically  in  the  left 


C  (m/s)  Distance  (m)  x  lO  ‘ 

Fig.  2  Result  of  the  ray  simulation.  The  average  vertical  profile  of  sound  speed,  used  in  the 
simulation  is  calculated  from  all  the  CTD  data  and  shown  at  the  left  of  the  figure.  The  bottom 
reflected  rays  with  reflection  numbers  less  than  two  are  drawn.  At  the  lower  middle  of  the  figure,  the 
travel  times  are  plotted  against  the  launch  angles  at  T2.  The  ray  paths  with  travel  times  around 
23.9s,  24.0s  and  24.1s  are  named  Ray-1,  Ray-2  and  Ray-3,  respectively.  Ray-1,  Ray-2  and  Ray-3  are 
drawn  with  the  thin  solid  line,  the  thick  dashed  line  and  the  thick  solid  line,  respectively. 


panel  of  Fig.  3,  there  are  several  arrival  peaks  in  each  of  correlation  waveforms,  acquired  every  eight 
hours.  We  here  focus  on  three  arrival  peaks  distributing  in  the  travel-time  range  from  23.9s  to  24.1s. 
These  arrival  peaks  are  identified  as  ray  groups  for  the  first,  second  and  third  arrivals.  In  comparison 
with  Figure  2,  these  arrival  rays  are  likely  to  correspond  to  Ray-1 ,  Ray-2  and  Ray-3.  The  significant 
arrival  peaks  sometimes  exist  at  24.2s  and  more,  but  these  data  corresponding  to  rays  with  multi 
reflections  at  the  surface  and  being  less  coherent  are  not  used  in  the  following  analyses. 

All  the  correlation  waveforms  obtained  originally  at  the  stations  T1  and  T2  are  stacked  with  the 
transmission  time  going  upward  as  shown  at  the  right  panel  of  Fig.  3. 
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Fig.  3  Typical  correlation  waveforms  obtained  at  T1  and  T2  (left  panel)  and  the  stack  diagrams  for 
the  T1  and  T2  data,  re-sampled  every  three  days  (right  panel).  The  arrival  peaks  corresponding  to 
Ray-1,  Ray-2  and  Ray-3  are  pointed  out  with  the  downward  arrows. 

This  is  the  result  before  cloek  eorrection.  Arrival  peaks  are  visible  as  a  wave  paeket  on  the  eorrelation 
waveform.  The  wave  packets  are  shifted  continuously  with  the  transmission  time.  The  internal  clock  of 
T1  and  T2  makes  a  reverse  tendency  of  clock  drift:  the  T1  clock  goes  ahead  and  T2  clock  is  delayed. 
The  necessity  of  clock  correction  is  well  understood  as  seen  by  the  large  drift  of  the  internal  clock. 


Velocity  data  and  comparison  with  other  data 


The  correction  rate  of  internal  clock  are  determined  by  putting  the  travel  time  difference  for  Ray-1  into 
zero.  The  corrected  travel  time  differences  (At)  for  Ray-2  and  Ray-3  are  converted  into  the  range- 

averaged  velocities (m„ ) ,  obeying  the  following  formula  (Zheng  et  ah,  1997): 


u 


m 


2L 
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where  L  denotes  the  horizontal  distanee  between  the  stations  T1  and  T2.  The  range-averaged  veloeities 
for  Ray-2  and  Ray-3,  smoothed  through  the  14-day  running  mean  are  shown  with  the  time  plots  in  Fig. 
4.  During  the  beginning  of  May  to  the  middle  of  July,  the  average  veloeities  are  direeted  toward  T1 
from  T2  (positive)  and  vary  in  the  range  of  0-10  em/s  and  5-15  em/s  for  Ray-2  and  Ray-3, 
respeetively.  After  that,  they  start  to  deerease  suddenly  and  make  a  reverse  eurrent  (direeted  from  T1  to 
T2),  eentered  at  the  beginning  of  August.  During  the  reverse  eurrent.  both  the  veloeities  are  almost  the 
same. 

The  ADCP  data,  obtained  at  the  stations  T1  and  T2,  are  also  shown  in  Figure  4  with  the  time  plots.  The 
three  layered  veloeities  of  100-200m,  200-300m  and  300-400m  are  shown  for  the  T2-ADCP  data  while 
only  the  100-200m  layer  veloeity  is  shown  for  the  T1  -ADCP  data.  Both  the  ADCP  data  for  the  100- 
200m  layer  are  in  good  agreement.  The  T2-ADCP  data  shows  a  signifieant  vertieal  shear  of  eurrent 
until  the  middle  of  July.  However  this  large  shear  is  suddenly  diminished  in  the  period  of  reverse 
eurrent,  eentered  at  the  beginning  of  August  and  all  the  OAT  and  ADCP  veloeities  are  eoineident. 

Summary  and  discussion 

The  oeean  aeoustie  tomography  experiment  is  earried  out  at  the  northern  part  of  the  Luzon  Strait  with 
an  aim  of  the  long-term  measurement  of  oeean  eurrent.  The  timing  module  eomposed  of  GPS,  Quartz 
and  Rubidium  eloeks  is  well  operated^  but  the  expeeted  performanee  is  not  attained.  The  aeeurate 
analysis  of  the  long-term  eurrent  data  is  aeeomplished  under  the  new  eloek  eorreetion  method  in  whieh 
the  eorreetion  rate  is  determined  by  putting  the  travel  time  differenee  for  rays  passing  near  the  USC 
axis  into  zero.  Information  on  oeean  dynamies  that  eurrents  are  generally  very  weak  at  the  bottom  of 
the  permanent  thermoeline,  namely  the  axis  of  USC,  plays  an  essential  role  in  the 
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Date  of  2008 


Fig.  4  OAT  velocities  for  Ray-2  (thick  black  line)  and  Ray-3  (thick  red  line),  compared  with  the 
ADCP  velocities  at  100-2 00m  (thin  blue  line) ,  200-300m  (thin  green  line)  and  300-400m  (thin  pink 

line)  of  T2  and  at  100-200m  (thin  green  line)  at  Tl. 


realization  of  this  clock  correction  method.  Yuan  et  al.  (2008)  reported  the  result  of  the  moored  current 
meter  that  mean  currents  acquired  at  depth  800m  in  the  adjacent  of  the  tomography  site  are  as  small  as 
Icm/s. 

During  the  observational  period  from  May  to  August  2008,  the  Kuroshio  is  likely  to  be  located  at  the 
east  of  the  tomography  site,  and  then  the  OAT  measures  a  relatively  weak  current  less  than  30cm/s  at 
the  western  edge  of  the  Kuroshio.  On  the  other  hand,  the  present  OAT  experiment  is  performed  at  the 
final  phase  of  the  2007/08  La  Nina  (Pacific  Marine  Environmental  Laboratory  Web 
Site,  http://www.pmel.noaa.gov/tao/  jsdisplay/).  Qu  et  al.  (2004)  shows  from  the  high-resolution  ocean 
circulation  model  that  the  Luzon  Strait  transport  (LST)  becomes  maximum  during  the  El  Nino  year 
and  minimum  during 
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Fig.  5  Sequential  maps  of  the  satellite  altimeter  data,  taken  around  the  observation  site  on  (a)  July 
20,  (b)  July  30,  (c)  August  10  and  (d)  August  20.  The  transmission  line  T1  -  T2  is  drawn  with  a  solid 

square,  placed  at  the  south  of  Taiwan. 


the  La  Nina  year.  Furthermore  during  the  La  Nina  years  the  southward  shift  of  the  North  Equatorial 
Current  (NEC)  bifurcation  occurs  and  the  Kuroshio  is  less  intruded  into  the  Euzon  Strait  with 
increased  transport  (Qu  and  Eukas,  2003).  This  is  the  reason  why  the  tomography  site  is  located  further 
west  of  the  stream  axis  of  the  Kuroshio.  The  ocean  circulation  modeling  work  also  shows  that  the  EST 
takes  seasonal  variability  to  be  maximum  in  January-February  and  minimum  in  June- July  (Qu  et  ah, 
2000).  This  provides  an  additional  cause  that  the  Kuroshio  takes  the  eastward-shifted  stream  path 
during  the  acoustic  tomography  experiment,  done  in  late  spring  to  summer. 
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At  the  overall  view  during  the  beginning  of  May  to  the  middle  of  July,  the  OAT  velocities  averaged 
along  Ray-  and  Ray-  are  significantly  weaker  than  the  ADCP  velocities  except  the  first  half  of 
June.  This  is  caused  by  that  the  former  velocities  average  the  depth  range  from  the  200m  to  the  bottom 
(Ray-2)  and  from  the  surface  to  bottom  (Ray-3)  while  the  latter  velocities  average  the  depth  range 
from  100m  to  400m.  In  contrast,  the  OAT  velocity  produces  an  overestimate  of  about  5cm/s  during  the 
first  half  of  June.  During  one  month  from  mid  July  to  mid  August,  all  the  observed  velocities  are 
surprisingly  coincident,  making  the  prominent  reverse  flows  and  diminishing  the  vertical  shear.  It  is 
likely  that  a  cyclonic  (cold)  eddy  with  baroclinicity  exists  between  the  tomography  site  and  the 
Kuroshio  in  the  Philippine  Sea  to  generate  this  reverse  current.  Cold  eddies  on  the  Philippine  Sea  side 
of  the  Luzon  Strait  are  visible  in  the  sea  surface  height  anomaly  (SSHA)  data  on  July  20,  July  30  and 
August  10,  retrieved  from  the  Naval  Research  Laboratory  Web  Site 

(http://www7320.nrlssc.navy.mil/modas2d/)  (Fig.  5,).  Especially  the  cold  eddies  in  Figures  5  (a)  and 
(b)  make  a  barrier  to  the  intrusion  of  the  Kuroshio  into  the  South  China  Sea.  The  cold  eddies  are  also 
reported  to  be  found  in  the  surface  geostrophic  current  calculated  from  the  combined  data  of  the  mean 
dynamic  topography  and  the  SSHA  (Figure  11  of  D.  Yuan,  2006).  It  is  reported  that  the  mesoscale 
eddies  impinge  the  Kuroshio  east  of  Taiwan  at  an  interval  of  100  days  (Yang  et  ah,  1999).  The  cold 
eddy  may  force  the  Kuroshio  to  make  the  eastward  shifted  path  on  the  eastern  side  of  the  eddy  and  the 
reverse  current  on  its  western  side  after  the  collision  with  the  Kuroshio.  However  the  detailed 
dynamics  of  Kuroshio-eddy  interaction  requires  further  study. 

During  the  reverse  current  the  OAT  velocities,  which  include  the  current  from  the  surface  to  bottom, 
agree  well  with  the  ADCP  velocities  for  the  upper  400m.  This  means  that  the  cold  eddy  is 
characterized  with  a  baroclinic  structure  and  its  velocity  is  almost  diminished  at  depth  500m, 
corresponding  to  the  shallowest  point  of  Ray- 1 .  This  baroclinic  structure  of  cold  eddies  is  also  roughly 
supported  by  the  CTD  section  acquired  at  25N  (Zhu  et  ah,  2008). 
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